R2R - software to speed the depiction of aesthetic consensus RNA secondary structures

<p>Abstract</p> <p>Background</p> <p>With continuing identification of novel structured noncoding RNAs, there is an increasing need to create schematic diagrams showing the consensus features of these molecules. RNA structural diagrams are typically made either with general-purpose drawing programs like Adobe Illustrator, or with automated or interactive programs specific to RNA. Unfortunately, the use of applications like Illustrator is extremely time consuming, while existing RNA-specific programs produce figures that are useful, but usually not of the same aesthetic quality as those produced at great cost in Illustrator. Additionally, most existing RNA-specific applications are designed for drawing single RNA molecules, not consensus diagrams.</p> <p>Results</p> <p>We created R2R, a computer program that facilitates the generation of aesthetic and readable drawings of RNA consensus diagrams in a fraction of the time required with general-purpose drawing programs. Since the inference of a consensus RNA structure typically requires a multiple-sequence alignment, the R2R user annotates the alignment with commands directing the layout and annotation of the RNA. R2R creates SVG or PDF output that can be imported into Adobe Illustrator, Inkscape or CorelDRAW. R2R can be used to create consensus sequence and secondary structure models for novel RNA structures or to revise models when new representatives for known RNA classes become available. Although R2R does not currently have a graphical user interface, it has proven useful in our efforts to create 100 schematic models of distinct noncoding RNA classes.</p> <p>Conclusions</p> <p>R2R makes it possible to obtain high-quality drawings of the consensus sequence and structural models of many diverse RNA structures with a more practical amount of effort. R2R software is available at <url>http://breaker.research.yale.edu/R2R</url> and as an Additional file.</p

Integrating sequence and structural biology with DAS.

BACKGROUND: The Distributed Annotation System (DAS) is a network protocol for exchanging biological data. It is frequently used to share annotations of genomes and protein sequence. RESULTS: Here we present several extensions to the current DAS 1.5 protocol. These provide new commands to share alignments, three dimensional molecular structure data, add the possibility for registration and discovery of DAS servers, and provide a convention how to provide different types of data plots. We present examples of web sites and applications that use the new extensions. We operate a public registry of DAS sources, which now includes entries for more than 250 distinct sources. CONCLUSION: Our DAS extensions are essential for the management of the growing number of services and exchange of diverse biological data sets. In addition the extensions allow new types of applications to be developed and scientific questions to be addressed. The registry of DAS sources is available at http://www.dasregistry.org.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Variation in Symbiodinium ITS2 Sequence Assemblages among Coral Colonies

Endosymbiotic dinoflagellates in the genus Symbiodinium are fundamentally important to the biology of scleractinian corals, as well as to a variety of other marine organisms. The genus Symbiodinium is genetically and functionally diverse and the taxonomic nature of the union between Symbiodinium and corals is implicated as a key trait determining the environmental tolerance of the symbiosis. Surprisingly, the question of how Symbiodinium diversity partitions within a species across spatial scales of meters to kilometers has received little attention, but is important to understanding the intrinsic biological scope of a given coral population and adaptations to the local environment. Here we address this gap by describing the Symbiodinium ITS2 sequence assemblages recovered from colonies of the reef building coral Montipora capitata sampled across Kāne'ohe Bay, Hawai'i. A total of 52 corals were sampled in a nested design of Coral Colony(Site(Region)) reflecting spatial scales of meters to kilometers. A diversity of Symbiodinium ITS2 sequences was recovered with the majority of variance partitioning at the level of the Coral Colony. To confirm this result, the Symbiodinium ITS2 sequence diversity in six M. capitata colonies were analyzed in much greater depth with 35 to 55 clones per colony. The ITS2 sequences and quantitative composition recovered from these colonies varied significantly, indicating that each coral hosted a different assemblage of Symbiodinium. The diversity of Symbiodinium ITS2 sequence assemblages retrieved from individual colonies of M. capitata here highlights the problems inherent in interpreting multi-copy and intra-genomically variable molecular markers, and serves as a context for discussing the utility and biological relevance of assigning species names based on Symbiodinium ITS2 genotyping

Contribution of a Common Variant in the Promoter of the 1-α-Hydroxylase Gene (CYP27B1) to Fracture Risk in the Elderly

CYP27B1 encodes mitochondrial 1α-hydroxylase, which converts 25-hydroxyvitamin D to its active 1,25-dihydroxylated metabolite. We tested the hypothesis that common variants in the CYP27B1 promoter are associated with fracture risk. The study was designed as a population-based genetic association study, which involved 153 men and 596 women aged 65–101 years, who had been followed for 2.2 years (range 0.1–5.5) between 1999 and 2006. During the follow-up period, the incidence of fragility fractures was ascertained. Bone ultrasound attenuation (BUA) was measured in all individuals, as were serum 25-hydroxyvitamin D and PTH concentrations; 86% subjects had vitamin D insufficiency. Genotypes were determined for the –1260C>A (rs10877012) and +2838T>C (rs4646536) CYP27B1 polymorphisms. A reporter gene assay was used to assess functional expression of the –1260C>A CYP27B1 variants. The association between genotypes and fracture risk was analyzed by Cox’s proportional hazards model. We found that genotypic distribution of CYP27B1 –1260 and CYP27B1 +2838 polymorphisms was consistent with the Hardy-Weinberg equilibrium law. The two polymorphisms were in high linkage disequilibrium, with D′ = 0.96 and r2 = 0.94. Each C allele of the CYP27B1 –1260 polymorphism was associated with increased risk of fracture (hazard ratio = 1.34, 95% CI 1.03–1.73), after adjustment for age, sex, number of falls, and BUA. In transient transfection studies, a reporter gene downstream of the –1260(A)-containing promoter was more highly expressed than that containing the C allele. These data suggest that a common but functional variation within the CYP27B1 promoter gene is associated with fracture risk in the elderly

Trends in modeling Biomedical Complex Systems

In this paper we provide an introduction to the techniques for multi-scale complex biological systems, from the single bio-molecule to the cell, combining theoretical modeling, experiments, informatics tools and technologies suitable for biological and biomedical research, which are becoming increasingly multidisciplinary, multidimensional and information-driven. The most important concepts on mathematical modeling methodologies and statistical inference, bioinformatics and standards tools to investigate complex biomedical systems are discussed and the prominent literature useful to both the practitioner and the theoretician are presented

A Proposed Taxonomy of Anaerobic Fungi (Class Neocallimastigomycetes) Suitable for Large-Scale Sequence-Based Community Structure Analysis

Anaerobic fungi are key players in the breakdown of fibrous plant material in the rumen, but not much is known about the composition and stability of fungal communities in ruminants. We analyzed anaerobic fungi in 53 rumen samples from farmed sheep (4 different flocks), cattle, and deer feeding on a variety of diets. Denaturing gradient gel electrophoresis fingerprinting of the internal transcribed spacer 1 (ITS1) region of the rrn operon revealed a high diversity of anaerobic fungal phylotypes across all samples. Clone libraries of the ITS1 region were constructed from DNA from 11 rumen samples that had distinctly different fungal communities. A total of 417 new sequences were generated to expand the number and diversity of ITS1 sequences available. Major phylogenetic groups of anaerobic fungi in New Zealand ruminants belonged to the genera Piromyces, Neocallimastix, Caecomyces and Orpinomyces. In addition, sequences forming four novel clades were obtained, which may represent so far undetected genera or species of anaerobic fungi. We propose a revised phylogeny and pragmatic taxonomy for anaerobic fungi, which was tested and proved suitable for analysis of datasets stemming from high-throughput next-generation sequencing methods. Comparing our revised taxonomy to the taxonomic assignment of sequences deposited in the GenBank database, we believe that >29% of ITS1 sequences derived from anaerobic fungal isolates or clones are misnamed at the genus level

Internal Transcribed Spacer 2 (nu ITS2 rRNA) Sequence-Structure Phylogenetics: Towards an Automated Reconstruction of the Green Algal Tree of Life

L). Some have advocated the use of the nuclear-encoded, internal transcribed spacer two (ITS2) as an alternative to the traditional chloroplast markers. However, the ITS2 is broadly perceived to be insufficiently conserved or to be confounded by introgression or biparental inheritance patterns, precluding its broad use in phylogenetic reconstruction or as a DNA barcode. A growing body of evidence has shown that simultaneous analysis of nucleotide data with secondary structure information can overcome at least some of the limitations of ITS2. The goal of this investigation was to assess the feasibility of an automated, sequence-structure approach for analysis of IT2 data from a large sampling of phylum Chlorophyta.Sequences and secondary structures from 591 chlorophycean, 741 trebouxiophycean and 938 ulvophycean algae, all obtained from the ITS2 Database, were aligned using a sequence structure-specific scoring matrix. Phylogenetic relationships were reconstructed by Profile Neighbor-Joining coupled with a sequence structure-specific, general time reversible substitution model. Results from analyses of the ITS2 data were robust at multiple nodes and showed considerable congruence with results from published phylogenetic analyses.Our observations on the power of automated, sequence-structure analyses of ITS2 to reconstruct phylum-level phylogenies of the green algae validate this approach to assessing diversity for large sets of chlorophytan taxa. Moreover, our results indicate that objections to the use of ITS2 for DNA barcoding should be weighed against the utility of an automated, data analysis approach with demonstrated power to reconstruct evolutionary patterns for highly divergent lineages

From Structure Prediction to Genomic Screens for Novel Non-Coding RNAs

Non-coding RNAs (ncRNAs) are receiving more and more attention not only as an abundant class of genes, but also as regulatory structural elements (some located in mRNAs). A key feature of RNA function is its structure. Computational methods were developed early for folding and prediction of RNA structure with the aim of assisting in functional analysis. With the discovery of more and more ncRNAs, it has become clear that a large fraction of these are highly structured. Interestingly, a large part of the structure is comprised of regular Watson-Crick and GU wobble base pairs. This and the increased amount of available genomes have made it possible to employ structure-based methods for genomic screens. The field has moved from folding prediction of single sequences to computational screens for ncRNAs in genomic sequence using the RNA structure as the main characteristic feature. Whereas early methods focused on energy-directed folding of single sequences, comparative analysis based on structure preserving changes of base pairs has been efficient in improving accuracy, and today this constitutes a key component in genomic screens. Here, we cover the basic principles of RNA folding and touch upon some of the concepts in current methods that have been applied in genomic screens for de novo RNA structures in searches for novel ncRNA genes and regulatory RNA structure on mRNAs. We discuss the strengths and weaknesses of the different strategies and how they can complement each other

Barcoding a Quantified Food Web: Crypsis, Concepts, Ecology and Hypotheses

The efficient and effective monitoring of individuals and populations is critically dependent on correct species identification. While this point may seem obvious, identifying the majority of the more than 100 natural enemies involved in the spruce budworm (Choristoneura fumiferana – SBW) food web remains a non-trivial endeavor. Insect parasitoids play a major role in the processes governing the population dynamics of SBW throughout eastern North America. However, these species are at the leading edge of the taxonomic impediment and integrating standardized identification capacity into existing field programs would provide clear benefits. We asked to what extent DNA barcoding the SBW food web would alter our understanding of the diversity and connectence of the food web and the frequency of generalists vs. specialists in different forest habitats. We DNA barcoded over 10% of the insects collected from the SBW food web in three New Brunswick forest plots from 1983 to 1993. For 30% of these specimens, we amplified at least one additional nuclear region. When the nodes of the food web were estimated based on barcode divergences (using molecular operational taxonomic units (MOTU) or phylogenetic diversity (PD) – the food web became much more diverse and connectence was reduced. We tested one measure of food web structure (the “bird feeder effect”) and found no difference compared to the morphologically based predictions. Many, but not all, of the presumably polyphagous parasitoids now appear to be morphologically-cryptic host-specialists. To our knowledge, this project is the first to barcode a food web in which interactions have already been well-documented and described in space, time and abundance. It is poised to be a system in which field-based methods permit the identification capacity required by forestry scientists. Food web barcoding provided an effective tool for the accurate identification of all species involved in the cascading effects of future budworm outbreaks. Integrating standardized barcodes within food webs may ultimately change the face of community ecology. This will be most poignantly felt in food webs that have not yet been quantified. Here, more accurate and precise connections will be within the grasp of any researcher for the first time